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A non-intmsive velocity diagnostic based on laser induced fluorescence of the 5d 4 F 5/2 - 6p 4 D 5/2 singly 
ionized xenon transition zvas used to interrogate the exhaust of a 1.5 kW Stationary Plasma Thruster (SPT). 
A detailed map of plume velocity vectors was obtained using a simplified, cost-effective, non-intrusive, 
semiconductor laser based scheme. Circumferential velocities on the order of 250 m/s were measured which 
implied induced momentum torques of approximately 5xlO~ 2 N-cm. Axial and radial velocities were 

evaluated one mm downstream of the cathode at several locations across the width of the annular 
acceleration channel. Radial velocities varied linearly with radial distance. A maximum radial velocity 
of 7500 m/s was measured 8 mm from the center of the channel. Axial velocities as large as 16,500 m/s were 
measured. 

Nomenclature 

A fine structure constant accounting for nuclear 
magnetic moment, MHz 

B fine structure constant accounting for electric 
quadrupole interaction, MHz 
C fine structure constant, MHz 
c speed of light, m/s 

F total angular momentum quantum number 
1 nuclear spin quantum number 
/ total electronic angular momentum quantum 
number 

k Boltzmann's constant, J/K 
m mass of ion, kg 
T temperature, K 

v, ion velocity in direction of laser, m/s 
W F quantum state term energy including 
hyperfine splitting, MHz 
Wj quantum state term energy excluding 
hyperfine splitting, MHz 
Av Doppler shift, cm" 1 
Av d Doppler width (FWHM), cm" 1 
v 0 unshifted frequency of electronic transition, 
cm' 1 

K constant 

Introduction 

The suitability of ion accelerators with 
crossed electric and magnetic fields, such as the 
Russian stationary plasma thruster (SPT), for 
satisfying various propulsion mission 
requirements is currently being assessed. 

Spacecraft integration issues are a major 
consideration for such a system to be successfully 
employed on next-generation spacecraft. NASA 
has taken an active role in assessing integration 
issues pertinent to the SPT using a thruster 
provided by the Ballistic Missile Defense 


Organization (BMDO). 1 " 5 One specific issue 
relevant to spacecraft integration of a SPT is the 
possible momentum torques generated by 
circumferential ion motion. Sputtering of 
spacecraft surfaces by ion bombardment also needs 
to be evaluated. 

The direction and magnitude of ion velocity 
vectors can be used to perform this type of 
spacecraft integration assessment. Such data 
would also aid in a detailed performance 
assessment because the spatial, integration of the 
density weighted velocity field determines the 
thrust vector. For example, the measurement of 
velocity in directions other than the direction of 
the thrust permits an evaluation of thrust 
inefficiencies due to beam divergence. 
Additionally, a detailed understanding of the 
flow field in the region of critical engine surfaces 
can provide information on sputtering due to 
impingment. 

This paper describes the direct measurement 
of the ion velocity at the exit of a SPT based on a 
high-spatial-resolution, non-intrusive optical 
diagnostic. To date the use of non-intrusive 
optical diagnostics applied to SPT's consisted of 
emission spectroscopy measurements of excited 
state populations in the exit plane of the SPT- 
100. 5 In this investigation the optical technique 
employed to measure ion velocities was laser 
induced fluorescence (LIF) using a semiconductor 
laser. LIF is particularly well suited for 
measurements near the thruster exit where 
survivability of intrusive probes is an issue and 
where steep property gradients make high 
spatial resolution advantageous. Measurements 
of the velocity field at the thruster exit are 


desirable because possible facility interactions 
are minimized at this location. These results can 
then be extrapolated to determine the ion 
impacts at various distances from the thruster. 
The use of a semiconductor laser enables 
implementation of a such a diagnostic for a 
fraction of the cost of traditional systems. 

Theory 

Velocity measurement using laser induced 
fluorescence is based on the Doppler effect. Laser 
induced fluorescence can be described as two 
separate processes: photon absorption and 
subsequent spontaneous emission. Absorption 
occurs when the interaction between an incident 
photon and an ion results in the promotion of a 
bound electron from one energy quantum state to a 
higher energy quantum state. This may take 
place when the photon has the same energy as 
the electronic transition. The excited state is not 
stable and will spontaneously decay while 
emitting a photon in an arbitrary direction. 

If the ion is moving relative to the laser 
beam, the laser photon energy in the moving 
frame of reference is Doppler shifted by an 
amount proportional to the relative velocity 
between the ion and laser. The change in photon 
energy (or laser frequency) needed to be resonant 
with the electronic transition of a moving ion is 
referred to as a Doppler shift and can be related 
to the resonant frequency for a stationary ion, the 
component of the ion's velocity in the direction of 
the laser beam, and the speed of light according 
to: 

Av = v 0 • (v, /c ) (1) 

This velocity diagnostic concept has successfully 
been implemented by several investigators in 
high velocity plasma flows. 6 " 8 

In this investigation a non-resonant scheme 
was implemented (ie. spontaneous decay of the 
upper state to a state other than the lower laser 
coupled state). The optically allowed 5d 4 F 5 /2- 
6p 4 D 5/2 transition of singly ionized xenon which 
has an energy of 11976.72 cm" 1 (834.7 nm) was 
stimulated while monitoring spontaneous 
emission from the 6p 4 D 5 /2 - 6s 4 P3/ 2 transition at 
18447.95 cm" 1 (541.9 nm). The laser coupled 
transition at 11976.72 cm" 1 was chosen for 
convenience as it was easily accessible using a 
semiconductor laser. 

Additional information such as kinetic ion 
temperatures and magnetic field strengths may be 
determined from fluorescence excitation line 
shapes. This type of analysis requires a detailed 


line shape model which should include isotopic 
effects, hyperfine splitting, Zeeman splitting, 
saturation broadening due to excessive laser 
power, and other broadening mechanisms 
associated with the plasma such as Doppler 
broadening, resonance broadening, Stark 
broadening, lifetime broadening, and transit time 
broadening. For the conditions in the SPT plume, 
Doppler broadening should dominate, so a model 
incorporating isotopic effects, hyperfine 
splitting, and Doppler broadening was 
constructed. 

Naturally occurring xenon is composed of 
seven isotopes. For each of these isotopes the 
energy of the measured transition is slightly 
shifted. The magnitude of the isotopic shift can 
be determined based on the change in the center of 
gravity of the ion, the change in the nuclear 
radius of the ion, and the degree of interaction 
between the ion's electrons. The change in the 
center of gravity can accurately be calculated. 
Calculations of the other effects requires accurate 
wave functions presently not available for xenon. 
Because of this, experimentally determined 
isotope shift values were sought, though they 
were not available for the investigated 
transition. In the absence of these data, 
published data for a similar transition, the 5d 
4 D 7/ 2 - 6p 4 P5/2 transition at 16521.24 cm" 1 (605.1 
nm) were used. 


Table I: Isotope Shifts 


Isotope 

% natural 
abundance 9 

Isotope Shift 
(rel to Xe 132 ) MHz 10 

Xe 132 

26.89 

0 

Xe 129 

26.44 

162 

Xe 131 

21.18 

99 

Xe 134 

10.44 

80 

Xe 136 

8.87 

12 

Xe 130 

4.08 

-82 

Xe 128 

1.92 

-145 


Hyperfine splitting occurring in the odd 
numbered isotopes gives rise to additional 
structure of the excitation line shape. Xe 129 has a 
nuclear spin quantum number I of 1/2 which splits 
the upper and lower states into two states. This 
results in four transitions. Xe 131 has a nuclear spin 
quantum number of 3/2 which splits both the 
upper and lower states into four individual 
states. Transitions with AF = ±1,0 are allowed, 
except from F= 0 to 0. The total angular 
momentum quantum number F takes values from: 

F = /+/,7+/-l I /-/I (2) 


2 


where J is the total angular momentum of the 
electrons. This results in ten allowed transitions 
which are shown in Fig. 1. 
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Fig. 1 Hyperfine Levels and Transitions for Xe 129 
and Xe 131 


The energy of the individual hyperfine 
levels is determined by modifying the term 
energy for the individual odd isotopes with 
correction terms due to the magnetic moment of 
the nucleus and the electric quadrupole coupling 
between the nucleus and electrons: 


k[J(J+1)+F(F+1)-I(I+1)] 2 (2F+1) 

~ F(F+1) (5 ' 

if AF = 1 

- K(]+F+M)(]+F-1)(J-F+1+1)(J-F-I) 

I- p (6) 

and if AF = -1 

-K(]+F+I+2)(]+F-l+l)(]-F+l)(]-F-l-l) 

( F+l) (7) 

where jc is an unknown constant. 

The width of each of the 20 individual 
components of the transition manifold may have 
been affected by any of the broadening 
phenomena previously mentioned. Additionally, 
each state may have been Zeeman split into 2F+1 
states if magnetic field effects are substantial. 
For this model only Doppler broadening was 
considered. Doppler broadening is the result of 
the spread in velocity for the species being 
investigated. This spread in velocity classes is 
generally attributed to the random thermal 
motion of the ions. This gives rise to a Gaussian 
shape for each of the transitions with a half 
width of: 


W f =W,+AC/2+B 


3/4C(C+l)-l(l+l )](]+!) 
21 ( 21 - 1 ) 1 ( 2 ]-!) 


(3) 


hv D = v 0 -yj' 


8kTln2 


( 8 ) 


where W r is the term energy including fine 
structure and Wj is the term energy accounting for 
everything other than hyperfine splitting. A is a 
fine structure constant which corrects for the 
nuclear magnetic moment and B is a fine structure 
constant that accounts for the electric quadrupole 
interaction between the nucleus and electrons. 
The quadrupole interaction vanishes for Xe 129 
which has a spherical nuclear charge 
distribution. The constant C can be calculated 
based on quantum mechanical considerations. 

C=F(F+1)-1(I+1)-J(J+1) (4) 

Typically fine structure constants are 
determined experimentally. These constants 
were unavailable for the transition considered in 
this investigation. 

As a result of these isotopic shifts and 
hyperfine splitting the 5d 4 F 5/2 6p 4 D 5 /2 
transition is split into 20 individual transitions. 
The relative intensities of the contribution from 
each isotope is determined based on the relative 
abundance of each isotope. The relative intensity 
of each hyperfine component can be calculated 
based on quantum mechanical considerations. For 
this transition (/= 5/2 to ]= 5/2) if AF = 0 


Apparatus and Procedure 

The investigation of ion velocity in the 
exhaust of the SPT-100 was conducted in a 
vacuum chamber measuring 1.5 m in diameter by 5 
meters in length. The pumping for this chamber 
consisted of four 0.82 m oil diffusion pumps (each 
rated at 30,000 1/s), a lobe type blower, and two 
piston type roughing pumps. The SPT was 
mounted in a 60 cm diameter spool piece 
separated from the main chamber with a 60 cm 
isolation valve. The smaller diameter of the 
spool piece permitted the entire optical system to 
be located outside the vacuum facility. The 
locations of the spool piece relative to the tank, 
the thruster relative to the spool piece, and the 
diffusion pumps are shown in Fig. 2. 
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The SPT-100 used for this investigation was 
fabricated by Fakel Enterprises in Russia and was 
provided by BMDO. It was previously used in 
several experimental investigations and a more 
complete description is given elsewhere. 1 A 
photograph of this thruster is shown in Fig. 3. 
The outer diameter of the annular discharge 



Fig. 3 SPT-100 

chamber was 10 cm. The inner diameter was 5.6 
cm. The exit plane of the cathode /neutralizer 
assembly was 1 cm downstream of the thruster 
exit plane. The thruster was run on commercially 
available research grade xenon (99.9995% pure) 
at a nominal operating condition of 300 V 
discharge voltage, 4.5 A discharge current, and 5 
mg/s propellant flow. The total gas flow was 
internally divided between the anode and 
cathode. The flow system used to supply the gas 
to the thruster consisted of a needle valve and 
flow meter downstream of a pressure regulator, 
and a 250 liter compressed xenon cylinder. The 
flow meter was calibrated using a volume 
displacement method. The power supply used to 
runthe discharge was a 600 V, 25 A DC power 
supply operated with an external 10 Ohm ballast 
resistor. A 100 mF capacitor was installed across 
the discharge to minimize the impact of plasma 
oscillations on the power supply. Additional DC 
power supplies were used to run the heater 
cathode during start up and the keeper prior to 


initiation of the main discharge. A simplified 
schematic of the electrical configuration and the 
instrumentation used to monitor operation is 
shown in Fig. 4. In all cases the thruster was 
allowed to reach steady state conditions before 
spectroscopic measurements were taken. 



Fig. 4 SPT Electrical Schematic 

The experimental apparatus for the laser 
induced fluorescence measurements consisted of 
the laser system, associated positioning and 
focusing optics, and the fluorescence detection 
system. The laser used was a commercially 
available continuous wave, single-mode, 100 mW 
GaAlAs semiconductor laser with a nominal 
center wavelength of 845 nm and a line width of 
approximately 15 MHz (0.0005 cm''*). Frequency 
stability of the laser was optimized by using a 
very stable, low noise battery-powered current 
source. Laser temperature was actively regulated 
using a thermoelectric heater/ cooler. Frequency 
tuning was accomplished by varying the injection 
current. 

During a typical scan the laser frequency was 
monitored using an evacuated scanning Michelson 
interferometer with a resolution of 0.01 cm' 1 . The 
mode structure of the laser was evaluated with a 
2 Ghz free spectral range confocal etalon. A 
anamorphic prism pair was used to circularize 
the elliptical beam and a Faraday isolator was 
used to minimize optical feedback to the laser 
cavity. The output laser power was monitored 
using a calibrated power meter. 

The laser beam was directed to the 
appropriate optical access port on the spool piece 
using front surface mirrors. A 400 mm focal length 
lens was used to focus the laser to a sub- 
millimeter diameter at the probe volume. Non- 
resonant fluorescence at 541.9 nm was collected 
using a 400 mm focal length lens aligned 
orthogonally to the laser beam. A fixed 
wavelength 0.5 m Czemy-Tumer monochromator 
was used as a spectral filter. The fluorescence 
was converted to a signal by a bialkali type 
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photo multiplier. The current produced by the 
photo multiplier tube was monitored as a voltage 
across a 1 kD load resistor. This voltage was 
measured using a digital lock-in-amplifier. The 
lock-in amplifier permitted phase sensitive 
detection of the fluorescence signal at 1.2 kHz, 
the frequency the laser beam was modulated at 
using a mechanical chopper. Phase sensitive 
detection was used to discriminate spontaneous 
emission at 541.9 nm from the non-resonant 
fluorescence signal at the same frequency. The 
entire laser induced fluorescence apparatus is 
shown schematically in Fig. 5. 



Fig. 5 Optical Schematic (for circumferential 
velocity determination) 

A series of three different tests were conducted to 
investigate ion velocities in the axial, radial, 
and circumferential directions. These different 
velocity components were measured by bringing 
the laser into the probe volume in three different 
directions. The circumferential velocity was 
measured at the exit of the discharge chamber in 
the center of the annular channel with the beam 
tangent to the annulus as depicted in Fig. 4. The 
radial and axial velocity components were 
measured at various positions at a location 11 mm 
downstream of the discharge chamber exit plane. 
The thruster was moved laterally with a 
computer-controlled, stepper-motor-driven, 
linear actuator system to accomplish this. The 
locations of these measurements are shown in Fig. 
6. The radial velocity component was measured 
by directing the laser in the minus r direction 
(also shown in Fig. 6). The axial velocity' 
component was determined in conjunction with 


radial velocity by bringing the laser into the 
probe volume at an angle of 60 degrees relative to 



Fig. 6 Measurement Locations for Axial and 
Radial Velocity Determination 

The laser was scanned over the fluorescence 
excitation spectrum at a rate of 0.002 Hz using a 3 
s lock-in-amplifier time constant. Several 
additional scans were taken at a scan rate of 
0.0002 Hz using a 30 s time constant. These longer 
scans improved the signal to noise substantially 
and were performed when high fidelity line 
shape data were desired. Some of the scans were 
repeated. 

Results and Discussion 

Circumferential Velocity 

The circumferential velocity component was 
the first to be investigated. Two fluorescence 
excitation spectra, measured with the laser 
propagating in opposite directions with respect to 
the plume, are shown in Fig. 7. The laser was 
oriented such that it was collinear to the 
azimuthal velocity component. 

Because the two measurements were taken at 
the same spatial location, in the center of the 
annular channel at the exit of the discharge 
chamber, the shift in velocity is twice the shift 
relative to a stationary source. The frequency 
centered between these two scans is 11976.72 cm' 1 . 
Energy levels published by Moore 11 obtained from 
emission spectra indicate an unshifted frequency 
of 11976.64 cm-1. The experimentally determined 
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unshifted line position was used through out this 
investigation as the unshifted frequency. 



Fig. 7 Fluorescence Excitation Spectra Showing 
Circumferential Velocity Shift 

The measured shift corresponds to a 
circumferential velocity of approximately 250 
m/s in the clockwise direction (looking towards 
the thruster). This direction of ion motion is 
consistent with the outward direction of the 
radial magnetic field. There was a 5xl0' 3 N-cm 
torque produced by the thruster assuming the 
thruster efflux was entirely composed of singly 
ionized xenon moving with this circumferential 
velocity. This corresponds to an equivalent force 
acting at a distance of 3.9 cm (the distance from 
the center of the thruster to the center of the 
annular discharge chamber) of approximately 2% 
of the measured thrust. 

This estimate should serve as an upper bound 
because the plasma is not entirely composed of 
singly ionized xenon and the radial magnetic 
field is likely near its peak value at the 
measurement location. However, even this 
simplified assessment of torque is consistent with 
values published by Kozubsky 12 who indicated 
the equivalent force was between 0.5% and 2% of 
the thrust. Other investigators have published 
higher estimates. Barnhart 13 indicates a lxlO' 2 
N-cm torque, or 3% of the thrust andVaughn 14 
suggests a torque as high as 8% of the thrust. 
These values suggest circumferential velocities 
somewhat larger than measured for singly 
ionized xenon. 

Radial Velocity 

The radial velocities shown in Fig. 8 were 
measured 11 mm downstream of the discharge 
chamber exit plane ( 1 mm downstream of the 


cathode). Each point was determined based on a 
curve fit of the measured fluorescence signal 
frequency variation. The 0 mm radial position 
denoted in Fig. 6 corresponds to the center of the 
annular channel. Coordinates in this case are 
positive going outward from the geometric center 
of the thruster (refer to Fig. 6). The variation of 
velocity with position indicates the necessity of 
high spatial resolution for accurate 
determination of radial velocity. 



Radial Position, mm 

Fig. 8 Radial Velocity Measured 11 mm 
Downstream of the Discharge Chamber Exit 
Plane. 

The flow was divergent with a near linear 
variation in radial velocity with radial 
distance. The position of zero radial velocity 
was determined to be several millimeters inward 
from the center of the channel. A radial velocity 
of 7500 m/ s was measured at a location 8 mm from 
the center of the channel. This was the largest 
radial velocity measured. Measurements were 
not made at larger radial distances, but the 
linear variation suggests that larger values of 
radial velocity may exist at larger radial 
distances. The magnitude of these radial 
velocities indicate that there are performance 
losses due to divergence in the exhaust for the 
SPT-100 operating at the single operating point 
tested. 
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Axial Velocity 

Axial velocities were determined in 
conjunction with the radial velocity at the same 
measurement positions. The Doppler shift was 
measured with the laser sixty degrees from the 
axial direction in the same plane as the radial 
measurements. Because the laser was not 
coincident with axial ion velocity, the shift was 
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due to the axial and radial velocities. The axial 
velocity contribution to the shift was determined 
by removing the radial velocity's contribution. 
Fewer axial velocity measurements were taken 
due to time constraints. These data are shown in 
Fig. 9. 



Radial Position, mm 

Fig. 9 Axial Velocity Measured 11 mm 
Downstream of the Discharge Chamber Exit 
Plane. 

The magnitude of the axial velocities were 
consistent with a mass average velocity of 
approximately 15000 m/s determined based on 
the measured specific impulse. The variation in 
axial velocity with discharge voltage was 
investigated by measuring the Doppler shift at 
the center of the channel while the thruster was 
operating at 250, 300, and 350V. Assuming there 
was no change in the direction of the velocity 
vector with discharge voltage the axial velocity 
increased from 12700 m/s at a discharge voltage 
of 250 V to 15600 m/s at 300 V to 15800 m/s at 350 
V. These values of the velocity varied less with 
the discharge voltage than the measured thrust 
reported by Brophy et. al. 15 

The axial and radial velocity measurements 
were used to construct velocity vectors. These 
vectors are depicted in Fig. 10. From this, the 
divergence of high energy ions was apparent. 
Data of this type may be of use for evaluating 
integration of such thrusters onto spacecraft 
because these high energy ions could cause 
significant sputter damage if they were to come in 
contact with a surface. Additional data of this 
type may allow a determination of zones of 
influence without relying solely on costly, time 
consuming direct impingement tests. 



Fig. 10 Velocity Vectors 11 mm Downstream of 
the Discharge Chamber Exit Plane. 


Plasma Parameters 

The shape of the fluorescence excitation 
spectrum was evaluated to assess additional 
plasma parameters. Saturation behavior of the 
fluoresence transition was investigated and 
several non-saturated line shape measurements 
were recorded. Fig. 11 shows a measured line 
shape and a curve fit based on the previously 
discussed model. The scan was taken during an 
axial velocity determination. Each of the 20 
individual transitions included in this curvefit 
were Doppler broadened an amount corresponding 
to 20,000 K. Published fine structure constants 
were used for the 6p state. 10 Values for 5d state 
where unavailable and were estimated based on 
the measured line shape. The temperature 
corresponding to the line shape measured in the 
azimuthal velocity determination was on the 
order of 800 K. This is easily seen by comparing 
the scan widths shown in Figs. 7 and 11. 



Laser Frequency, cm- 1 

Fig. 11 Fluorescence Excitation Spectrum 
Measured During an Axial Velocity 
Determination and Analytic Fit. 
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The temperature used in the model to descibe 
the line width reflects a velocity spread in the 
laser direction. The experimental data indicate 
that the spread of velocity depended on 
direction. The spread in circumferential velocity 
was attributed to the ion’s thermal motion 
because no circumferential potential gradients 
were expected. Therefore, the 800 K temperature 
used to model the scan taken in this direction was 
interpreted as the ion's kinetic temperature, a 
temperature near that of the discharge chamber 
walls. The comparatively large temperature 
needed to model the axial line shape indicated a 
spread in velocity that may be attributed to a 
variation in accelerating voltage for the ions. If 
the exhausting ions were created in a region of 
space with a potential difference there would be 
a spread in axial velocity corresponding to this 
potential difference. The velocity spread 
indicated by the scan shown in Fig. 11 corresponds 
to an approximately 3.4 eV variation in ion 
energy (twice the 20000 K measured at 60 degrees 
from the axial direction). The largest velocity 
determined irrespective of direction was nearly 
18000 m/s corresponding to an ion energy of 
approximately 220 eV which is large with 
regard to an 3.4 eV energy spread. 

Conclusions 

A non-intrusive, laser induced fluorescence 
diagnostic was used to measure the velocity field 
for singly charged xenon in the exhaust of a SPT. 
Axial, radial, and circumferential velocities 
were measured and velocity vectors w r ere 
reconstructed from these data. These data 
demonstrate that detailed information with 
high spatial resolution can be obtained which 
can be used for performance assessment and 
spacecraft integration of these thrusters. 

Momentum torques caused by ion trajectories 
with a circumferential velocity component were 
estimated. Based on measured circumferential 
velocities, this torque was conservatively 
estimated to be 5xl0‘ 2 N-cm. Radial velocity 
measurements indicate a relatively divergent 
flow field. Axial velocities consistent with the 
measured thrust were determined. Ion velocity 
vectors in the exhaust of the SPT were 
reconstructed from the measurements of the 
individual velocity components. The magnitude 
of the ion velocity vector indicates that ion 
energies as high as 220 eV are present in the 
exhaust with an approximately 3.5 eV energy 
spread. This spread in energy was determined 
based on the measured fluorescence excitation 


spectrum and an analytic model constructed to 

simulate this transition. 
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